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Abstract 


k wave number, Eq. (35) 


A finite element model was developed to solve 
for the acoustic pressure field In a nonhomogeneous 
region. The derivations from the governing equa- 
tions assumed that the material properties could 
vary with position resulting In a nonhomogeneous 
variable property two-dimensional wave equation. 

This eliminated the necessity of finding the bound- 
ary conditions between the different materials. 

For a two media region consisting of part air (In 
the duct) and part bulk absorber (In the wall), a 
model was used to describe the bulk absorber proper- 
ties In two directions. An experiment to verify 
the numerical theory was conducted In a rectangular 
duct with no flow and absorbing material mounted on 
one wall. Changes In the sound field, consisting 
of planar waves was measured on the wall opposite 
the absorbing material. As a function of distance 
along the duct, fairly good agreement was found In 
the standing wave pattern upstream of the absorber 
and In the decay of pressure level opposite 
the absorber. 

Nomenclature 

A amplitude 


k xn axial modal wave number, Eq. (33) 

L dimensionless length, L'/b a 

L t,n non-dimensional parameter, Eq. (6) 

L v n L v p non-dimensional parameters, Eqs. (17 
* * and (18) 

m mode number, Eq. (38) 

N m number of modes 

P dimensionless pressure, 

P'fx.y.tJ/piactf 

Pr Prandtl number 

p dimensionless pressure, P( x,y ,t)/e iwt 

R dimensionless gas constant, R'To a /Col 

T n experimental constant 4 

T 0 dimensionless temperature, Tq/To 3 


A n 

A n 

BS 

b' 


mode amplitude of plus going entrance t 
waves, Eq. (32) 

u 

mode amplitude of reflected negative 
going entrance waves, Eq. (32) 

v 

mode amplitude of plus going exit waves, 

Eq. (36) 

x 

characteristic duct height, chosen 

bi y 


dimensionless 

dimensionless 
u ' /c oa 

dimensionless 

v’/ci a 

dimensionless 

dimensionless 


time, Eq. (1) 

axial acoustic velocity, 

transverse acoustic velocity, 

axial distance, Eq. (1) 
transverse distance, Eq. (1) 


b a dimensionless entrance height, 

b>' 

bfo dimensionless exit height b^/h' 

c e effective speed of sound, Eq. (9) 

c 0 dimensionless speed of sound, -J yRT 0 , 

Eq. (5), (Co/ C ; a = 1) 

Co speed of sound, ^ yR ' Tq 

d fiber diameter, d'/b a 

E property tensor, Eq. (30) 

f l frequency 

K dimensionless heat transfer parameter, 

K'ba/poa c oa» Eq. (20)* 


y specific heat ratio 

c property constant, Eqs. (25) and (26) 

C porosity 

v property constant, Eq. (27) 

t> 1 fluid kinematic viscosity 

p dimensionless density, Eq. (2)(p 0 = 1) 

a dimensionless viscous loss coefficient, 

• . • * i i 

a D a'Poa c oa 

a‘ viscous less coefficient, Eq. (13) 4 

g> dimensionless frequency, w'b a /Co a 

o' angular velocity 
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Subscripts 

a Inlet duct condition 

b exit duct condition 

1 1 th nodal point 

o ambient conditions 

x,y scalor vector components 

Superscripts 

' dimensional quantity 

approximate finite element solution 
Introduction 


The relatively high fuel economy available 
from propeller-driven aircraft has renewed Interest 
In high speed, highly loaded, multiple blade turbo- 
prop propulsion systems. The acoustic testing 1 of 
a scaled multiple blade propeller has been carried 
out In the 9- by 15-Ft Low Speed Wind Tunnel at NASA 
Lewis Research Center. Under simulated flight con- 
ditions In the wind tunnel, It Is desirable to meas- 
ure both acoustic amplitude and directivity to help 
fully characterize the propeller noise source. This 
Is not possible If reflections from the wind tunnel 
walls are Interfering with the direct sound. Conse- 
quently, the NASA Lewis 9- by 15-Ft Low Speed Wind 
Tunnel was acoustically lined to reduce tunnel 
reverberations . 

At the same time as measurements were made In 
the wind tunnel, analysis was done to provide data 
on how the measurements were affected by various 
acoustic parameters. In Ref. 2, Baumelster and 
Eversman used finite element theory to calculate the 
acoustic field of a propeller In a soft walled wind 
tunnel and compared these results to the radiation 
patterns of the same propeller In free space. They 
found that complicated pressure fields exist In a 
lined wind tunnel and that only over limited por- 
tions of the tunnel can measurements approximating 
the free field be obtained. 

The hard Helmholtz resonator cavity and a bulk 
fibrous absorber are the most common materials 
employed to attenuate sound. A bulk absorber called 
Kevlar was chosen as the wall absorbing material In 
the 9- by 15- Ft NASA Lewis Wind Tunnel. Design 
and construction of the wall liner was performed at 
NASA Lewis and Is documented In Ref. 3. In contrast 
to a locally reacting liner, such as the Helmholtz 
resonator shown In Fig. 1(a), the bulk extended 
reaction liner permits wave propagation In the axial 
direction, as shown In Fig. 1(b). 

In Ref. 2, the absorbing characteristics of the 
wind tunnel walls were modeled by applying the clas- 
sical admittance boundary conditions at the tunnel 
walls. Consequently, wave propagation In the axial 
direction In the bulk liner was not considered. The 
present Investigation wll 1 use finite element theory 
to develop a model to conveniently handle wave prop- 
agation In bulk material as well as In any Inhomog- 
eneous medium. The absorbing characteristics of the 
bulk material used In this model rely on the theo- 
retical development presented by Hersh. 4 


First, the geometric model and the appropriate 
governing equations and boundary conditions will be 
presented. Next, the finite element procedure for 
solving the problem will be discussed. Then, a 
numerical solution and an analytical solution to a 
test problem will be compared. Finally, the results 
of a validation experiment will be discussed and 
compared to theory. 

Geometrical Model 


In the finite element model of a tunnel with 
bulk absorbing materials along the wall, the present 
paper will focus on the Interaction of propagating 
duct modes with the wall lining as shown In Fjg. 2. 
In the uniform, Infinitely long entrance and exit 
regions with perfectly hard wall, the exact solution 
of the governing differential equations can be 
easily written In terms of the duct modes. In the 
central region which Includes both the duct and the 
fibrous absorbing region, the finite element analy- 
sis Is employed to determine the ressure field In 
this nonuniform region. Pressure waves propagating 
down the duct are either reflected, transmitted or 
absorbed by the nonuniform segment of the duct. 

Pressure mode reflection at the Inlet to the 
absorbing region and transmission at the outlet of 
the absorbing region are modeled by matching the 
finite element solution In the Interior of the cen- 
tral region to the known analytical eigen function 
expansions In the uniform Inlet and outlet ducts. 
This permits a multimodal representation accounting 
for reflection and mode conversion by the nonuniform 
absorbing section. This approach has been found to 
accurately model reflection and transmission 
coefficients . 5 


Governing Equations 

The governing equations are the state, continu- 
ity, and momentum linearized gas dynamics equations 
In the absence of a mean flow. In the fibrous mate- 
rial, the Hersh form 4 of the governing equations 
will be employed. By treating the bulk material as 
a momentum and thermal sink, Hersh modified the con- 
ventional linearized gas equation such that new 
complex propagation constants would account for 
acoustic energy absorption In the material. The 
development was semi-empirical since an adjustable 
constant was employed; however, the dimensional 
parameters developed appear to accurately follow 
experimental trends. Different sets of constants 
were determined for the fibers oriented perpendicu- 
lar to or parallel to the Incident sound particle 
velocity field. Consequently, some properties In 
the following equation will be subscripted with 
either x or y to Indicate anisotropic behavior. 

The Hersh model 4 was developed for acoustic 
propagation In one dimension. Herein, the model 
will be extended to two-dimensional propagation. 


The differential equations used herein are put 
In non-dimensional form form using 
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In the foregoing equations, the prime, , Is used to 
denote a dimensional quantity and the umprlmed 
defines a dimensionless quantity. This convention 
will be used throughout this paper. These and all 
other symbols used In the report are defined In the 
Nomenclature. Generally, conventional acoustic 
notation will be used so that the significance of 
most symbols are readily recognized. 



Thus, the equation of state can now be written as 
P(x.y.t) = cg P (x,y,t) (10) 


Equation of State 


Continuity Equation 


The equation relating acoustic pressure and 
density used In the bulk absorber was derived In 
Ref. 4 using a one-dimensional model for the heat 
•transfer between the sound field and the fibers. 

The heat transfer equation was used to determine a 
relationship between temperature and density per- 
turbations In the fluid surrounding the fibers. 

This relationship and the linearized perfect gas 
law yielded the following relationship between pres- 
sure and density perturbations which are assumed to 
hold In two-dimensions: 


P(x,y,t) 



(3) 


where the dimensionless frequency is defined as 

oj'ba' 

w = — ( 4 ) 


oa 


The parameters K Is defined In Ref. 4 as 

JL_ , ir ( _L 

p o w P r\w‘d ,2 j^ C J L t.n 


where 


In 


|Eq. (21) Ref. 4} (5) 

T 


5.4^1 - CH + 3.94(1 - C) 3 ] 

{Eq. (18b) Ref. 4| (6) 


In Eqs. (3) and (5), the symbol K denotes a 
heat transfer parameter that Is connected to the 
properties of the fluid medium and the fiber 
material (Eq. (20)). 4 


The two limiting cases of the state equation 
correspond to the Isothermal and adiabatic relation- 
ships between the pressure P and density p. At 
very low frequencies wherein K/p 0 o> » 1, 

P(x,y,t) « RT 0 p(x.y.t) (7) 

which represents the Isothermal relationship. At 
very high frequencies, wherein K/p 0 u « 1, 


P(x,y,t) * YRT 0 p(x,y,t) (8) 


the adiabatic relationship Is recovered. Further 
Interpretation of Eq. (3) can be found In Ref. 4. 

For simplicity. In the general case as given by 
Eq. (3), the effective speed of sound Is now defined 
as 


The continuity equation In the bulk material 
can be written as 


ap(x.y) p o 3U 3V _ n . 

at c ax c ay u ' 

where £ represents the porosity of the medium. In 
the main Inlet channel without fibers, £ has a 
value of unity. In the fibrous material £ Is a 
fraction less than one whose value depends on the 
packing density of the material. In this continuity 
equation, the particle mass flow is conserved where 
the air (or any gas) density p 0 Is constant In the 
medium and In the air gaps surrounding the fibers. 

Substituting the equation of state, Eq. (10), 
Into Eq. (11) and taking the derivative with respect 
to time yields 


_a£u 

atax 





Momentum Equations 


( 12 ) 


The x and y momentum equation can be 
written as 


p o au 
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(14) 


Here, the momentum equations have been empirically 
modified with a viscous loss coefficient a 
(Eq. 13, Ref. 4) resulting from skin friction, 
shear stresses and drag. 
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where V n and V p are arbitrary constants 
based on measurements and the parameters t Vfn 
and L VfP are those used In Ref. 4 


and 



] 

16^1 - C[l * 14.75(1 - C) 3 ] 
J Eq. 12, Ref. 4 1 


(17) 


d 


1 

3.94(1 - C)°' 413 [l + 27(1 - C) 3 ] 

{Eq. 11. Ref. 4j (18) 
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This Is an extension of the one-dimensional model 
used In Ref. 4 which used a viscous loss coeffi- 
cient corresponding to 


= <? x + o y j Eq. 13, Ref. 4| 


(19) 


The x momentum Eq. (13) can be rewritten In the 
form 


!p au 
C at 


ap _ _ j(_ au 
ax clu> at 


( 20 ) 


where u has been assumed to be a harmonic function 
of time. Thus, 


au 


-1 


dP 

^ax 


3t ~(t o . fx V 

\C cw 

Taking the x derivative of Eq. (21) yields 
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afu_ 

axat 


a_ 

ax 


Wave Equation 


Lc ci»J 


ap 

ax 


( 21 ) 


(22) 


Substituting Eq. (22) and the similar y form 
of the momentum equation Into the modified form of 
the continuity Eq. (12) yields 
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(23) 


Assuming, the pressure perturbation P to be 
harmonic In time, 

P(x,y,t) = p(x,y) e 1fc,t 
Eq. (23) becomes 
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For simplicity, let 


p 0 

10 A _X 

c x = c ' 

P 0 

, = -a . i -Jt 

y c o 


p o c e 


Thus, the wave equation becomes 
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(25) 

(26) 

(27) 

(28) 


V • j E • Vp)» pu 2 p = 0 (29) 

where the property tensor E Is represented by 

'M 

E =1 (30) 

0 M 

£ y/ 

and the vector product of the tensor E and a 
vector vp follows the common definition 
(Eq. A. 4-19, Ref. 6) 

Uniform Duct Analytical Solutions 

The analytical solutions of Eq. (28) for wave 
propagation In a uniform hard wall duct having an 
anechoic entrance and exit will be employed to give 
the termination boundary condition for the finite 
element region. The analytical solution for pres- 
sure waves traveling between parallel hard plates 
where the boundary condition Is 


|^=0 at y =0 and y » 

Is given as (see for example Ref. 7, p. 504) 
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time dependance used here, the 


(32) 


For the e 
-Ik x 

A + e xn term represents a wave propagating 

^ 1 k x 

- xn 

In the positive x direction while the A n e 

represents a wave moving In the negative x 
direction. 

The axial wave number k* n In Eq. (23) Is 
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(33) 


(34) 


(35) 


Equation (28) represents the governing wave equation 
to be solved by finite element theory. 

For later use In applying the finite element 
theory, It Is convenient to express Eq. (28) In 
vector form, 


The modal expression represented by Eq. (32) 
has been truncated to a total of N m modes of the 
Infinite number possible. Thus, a total of N m 
unknown modal amplitudes Aj, A?, — Apj m have 
been Introduced. The A + terms will be assumed as 
known. N m constraint equations will be required 
to determine each of these unknown reflection coef- 
ficients. The equations used to define these coef- 
ficients will be Introduced In the following section 
on boundary conditions. 

A similar solution exists at the exit, except 
only positive going waves are considered 


4 



Velocity Continuity 


P b (x,y) -]£ B n + cos y) e 1k><nX ( 36 > 

where k xn In Eq. (36) Is based on characteristic 
duct height of exit b^. 

Boundary Conditions 

A variety of boundary conditions will be used 
In the finite element solution of Eq. (28) for the 
model problem which Is displayed In schematic form 
In Fig. 4. Each of the required conditions will now 
be briefly discussed. 

I nput Condition 

The analysis assumes a given number N m of 
propagating Aj!j modes (see Eq. (32)). These modes 
effectively set the level of the scalar pressure 
field In the finite element region and can be viewed 
as the equivalent Dlrlchlet boundary conditions 
required for elliptic boundary value problems as 
defined by Eq. (28). 

Pressure Continuity at Inlet and Outlet 


In addition to pressure continuity, the axial 
acoustic velocity must also be continuous across the 
Interface to the finite element region. 


u = u (x = 0; 0 < y < b ) (40) 

a a 

Again, u a Is the modal representation of the 
scalar axial velocity field In the analytical solu- 
tion and u represents the approximate finite ele- 
ment solution. Using the momentum equations 
(Eq. (21)) to express the axial velocity In term of 
the pressure field yields 


1_ -^A L_ 

p a 3X = e x ax 


(41) 


where again an e +1wt time dependance has been 
assumed. The weak form of the finite element solu- 
tion will be employed In the solution of this prob- 
lem. In this form, a contour Integral term will be 
developed which will contain a natural boundary con- 
dition of the form vp • n where n represents the 
unit outward normal. Equation (41) can now be 
generalized to 


In general, the scalar pressure field Is con- 
tinuous across an Interface except where sources 
are present. Thus, the boundary condition at the 
Interface between the entrance duct and the finite 
element region requires 

p a - p (x , 0; 0 < y < b a ) (37) 

where p a Is the modal representation of the 
scalar pressure field In the analytical Inlet 
region given by Eq. (32) and p represents the 
finite element approximation for pressure at the 
Interface. The hat over p Implies an approximate 
finite element numerical solution to the true 
solution . 

At the Inlet Interface, shown by the dashed 
line In Fig. 4, the scalar pressure p a In the ana- 
lytical region given by Eq. (32) must match the pres- 
sure field defined by the finite element modal points 
along the boundary Interface. A weighted residual 
approach was used herein with the weighting function 
equal to the eigen functions. 

| 3 [P a (y) - P(y) ] C05 ( (m ~b a )lfy ) dy x 0 = 0 

(38) 

(N m equations m = 1,2,3, - - - N m ) 

Equation (38) represents N m separate equations; 
one for each coefficient defined In Eq. (32). The 
symbol m has been Introduced for the mode number 
equation to make It distinct from the multiple n 
mode number that make up the P a and analyt- 
ical functions. 

At the exit, 

b [p b (y) - p(y) 3 cos b b 1),rY ^ d y x _l = 0 

(39) 


-« x 3p a 

vp • n * — - - - (Inlet) (42) 

P a ax 

and at the exit 

+e ap. 

vp • n = — - — ^ (exit) (43) 

p 5 9x 

In general, the gradient of the pressure field Is 
continuous except where there Is a step change in 
density. The sign change In Eqs. (42) and (43) 
comes directly from the directional change of the 
unit outward normal . 

Hard Wall Boundary Condition 

At the hard walls shown In Fig. 4, the trans- 
verse acoustic velocity Is zero. Again, using the 
momentum equations to relate the acoustic velocity 
to the pressure fields requires 

vp • n = 0 (44) 

FINITE ELEMENT THEORY 

The finite element formulation of the Inhomog- 
eneous wave equation Is now generated by using the 
weighted residual approach with the Galerkln approx- 
imation to obtain an Integral form of the variable 
property wave equation over the whole (global) 
domain . 

The continuous domain D Is first divided Into 
a number of discrete areas as shown In Fig. 2. The 
particular pattern chosen has been found to give 
accurate results for a minimum number of elements. 8 

In the classical weighted residual manner, the 
pressure field Is curve fitted In terms of all the 
unknown modal values p^(x^,y^). The finite element 
aspects of converting Eq. (28) and the boundary con- 
ditions Into an appropriate set of global difference 
equations can be found In text books (e.g., Ref. 9) 
as well as Ref. 10 and for conciseness will not be 
presented herein. In Ref. 10, an exact finite ele- 
ment analog of Eq. (28) for electromagnatlc propaga- 
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1 1 on has been solved by Galerkin finite element 
theory with linear elements. 

Experimental Apparatus and Procedure 

A test apparatus was constructed to verify the 
numerical theories for a tunnel having a fibrous 
absorbing wall and no mean flow. The general acous- 
tic duct system described In Ref. 11 was modified 
Into the simple no-flow apparatus as shown schemat- 
ically In Fig. 3. The system was adapted to simu- 
late plane wave propagation In an Infinite duct. 

The following paragraphs describe the acoustic sys- 
tem including components and general operation. 

Source 

Sound was generated by a 120 W, 300 Hz to 6 kHz 
driver. In the experiment reported herein, a random 
noise generator provided a signal to the amplifier 
which In turn powers the acoustic driver. The duct 
temperature was at 28 °C. Since the first nonplane 
mode begins propagating at a frequency of 1701 Hz, 
the choice of 1560 Hz for the data frequency guaran- 
tees that only plane waves will propagate In the 
straight portion of the duct far from the fibrous 
absorption area. 


versed In the axial direction In Increments of 
0.953 cm (3/8 In.). The microphone holder plates 
were fabricated with an Insulator In each hole, so 
that the microphones were electrically Insulated 
from the test plate. During a measurement, the 
microphone's membrane was positioned flux with the 
surface of the test plate (protective shields were 
removed) . 

The response of these microphones was flat to 
25 kHz. a standard acoustic source was used to 
check the system amplification. The measured (nor- 
malized) r.m.s. acoustic pressures are shown In 
Table 1. The lengths are normalized to the duct 
height (3.81 cm) and the rms acoustic pressures are 
normalized by microphone (A). The flberous material 
begins at x/3.81 of 4.47 and ends at 8.167. 

Discussion of Results 

For theory and code validation, the finite 
element solution Is first applied to a case where 
an exact analytical solution exists. Next, the 
experimental results are compared to theoretical 
predictions . 

Reflection and Transmission With Normal Incidence 


Te st Section 

The sound travels through a 10 by 3.81 cm rec- 
tangular test section. This 1.4 m long section 
consists of 16 flat detachable plates (eight on top 
and eight on bottom). The fibrous absorbing test 
piece, to be described shortly, was Inserted In 
place of the bottom flat piece at location five. 

An exponential horn was also attached to the down- 
stream end of the test section. For these experi- 
ments, the horn was fitted with an acoustic foam 
wedge to absorb the Incident sound to approximate 
an anecholc termination. This anecholc termination 
simulates the Infinitely long exit duct assumed In 
the analysis. 

Fibrous Test Piece 

The fiber test piece was constructed from epoxy 
and metal. A 4.6 cm deep by 14.36 cm long and 
9.2 cm wide cavity was filled with Kevlar fibers 
with a overall density of 0.00792 g/cc, a mean fiber 
diameter of 0.61254 mm and with an effective poros- 
ity c of 0.9945. A four mesh metal screen covered 
the face as well as compressing the Kelvar In the 
cavl ty . 

Microphone Installation 


The first case considers a step change In mate- 
rial density from a non-dimensional value of 1 to a 
value of 4 at an axial position of x equal to 0.25 
Inside the finite element grid. In this case, the 
porosity Is unity, and the viscous loss coefficients 
<? x and <?y are assumed to be zero so that e x 
and e y take on the real values of 1 and 4, respec- 
tively. The y parameter was assumed to be 1 and 
the dimensionless Incident frequency was assumed to 
be 2ir . 

As shown In Fig. 5, the finite element and 
exact analytical theories (Ref. 12, p. 83) are In 
excellent agreement for the rms pressure. The 
reflections from the Interface between the two den- 
sity changes are clearly represented by the standing 
wave pattern ahead of the Interface. Observe that 
the magnitude of the pressure Increases Inside the 
material . 

Experimental Axial Pressure Profile - Hersh 
Properties 

The normalized rms pressure from Table 1 are 
plotted against axial distance In Fig. 6 and com- 
pared to the finite element solution of Eq. (28). 

The material properties were evaluated using Hersh's 
correlation 4 with 


Two 0.64 cm (1/4 In.) diameter condenser micro- 
phones were used to determine the acoustic field. 

The microphone designated (A) was used to monitor 
the source strength of the signal, while microphone 
(B) was used to determine the pressure and phase 
angle along the flat surface above the fibrous 
absorption test piece. The output of each micro- 
phone was analyzed to determine the r.m.s. pressure 
at a test frequency of 1560 Hz. The flat microphone 
test holder plate Is shown In Fig. 4 along with the 
fibrous absorption piece mounted In the duct. Three 
microphone holder plates were mounted above the 
fibrous absorption test piece at locations 4, 5, 
and 6 shown In Fig. 3. This aluminum microphone 
plate was fabricated with 0.64 cm (1/4 In.) diameter 
holes which allowed one microphone (B) to be tra- 


P r = 0.71 v 1 — 1 . 51 E-05 m 2 /sec 

•y : 1.4 d ' = 1 .25410-5 m 

Vp-1.0 V n = 0.44 T n - 1 .0 

and C - 0.9945 

and found to be 

( 44 ) 

(45) 

(46) 


c x = 1.0055 -10.0996 
Cy = 1 .0055 -10.1131 
v = 1.0004 -10.0479 


As seen In Fig. 6, the theory and experiment are 
In reasonable agreement. The magnitude of the 
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entrance standing wave Is In close agreement 
although a slight shift In the pattern Is seen. 

The standing wave pattern determines the reflected 
wave. The fall off of pressure In the central por- 
tion of the duct Is also handled reasonably well. 
However, the pressure rise at the exit Is over 
estimated although the proper trend Is predicted. 

The slight standing wave pattern In the exit 
duct data Is a result of reflections off the absorb- 
ing wedge In the test section horn. Since the 
theory assumes that there Is no reflected energy at 
the exit, a downstream standing wave pattern cannot, 
of course, arise In analytical predictions. The 
absence of reflected energy at the duct exit In the 
theory may account for the shift In the predicted 
standing wave pattern at the entrance. 

It can be shown that the standing wave pattern 
In a duct can be written as 



= 1 + R? + 2RrC0s(2kx + 0,) 
1 A * T 

(47) 

where Rf Is the magnitude of the termination 
reflection coefficient and Of Is the reflection 
coefficient phase angle. The standing wave pattern 
has a minimum value when 


c ° s(2lx .ln * V ' - 1 

(48) 

or 

2kx „,„ ‘ - (2 "- 1) * - 6 f 

(49) 


which applies when the axial coordinate gets more 
positive toward the termination. A measurement of 
the termination reflection coefficient at 1560 Hz 
shows that 6f = -0.45 ir, which results In a shift 
of the minima towards more positive values by about 

0.67 axial units. The data points agree with the 
direction of this shift but the measurements shown 
In Fig. 6 have about a 0.5 axial unit shift compared 
to the theory which Is only a 0.17 axial unit dif- 
ference. This difference Is within the accuracy of 
the microphone measurements. 

E xper i mental Axial Pressure Profiles - Modified 
H ersh Properties 

Other observations suggest that Hersh's empir- 
ical coefficient a could be underestimating vis- 
cous effects. Using an Impedance tube, Dahl and 
Rice 2 reported larger measured absorption than 
predicted with the Hersh modal In the 1500 Hz 
frequency range. Consequently, to check the sensit- 
ivity of the predictions to the viscous loss coef- 
ficient c, the Imaginary component of e has been 
Increased to 

e x = c y = 1 .0055 -10.21335 (50) 

As seen In Fig. 7, the theory and experiment are In 
better agreement at the exit. The present numerical 
theory may offer an attractive theoretical approach 
to evaluating viscous loss coefficients In bulk 
absorbers . 

Concluding Remarks 

The finite element model was developed to solve 
for the acoustic pressure field In a region that 
was nonhomogeneous . The derivations from the gov- 


erning equations assumed no mean flow and that the 
material properties could vary with position 
resulting In a nonhomogeneous variable property two- 
dimensional wave equation. This eliminated the 
necessity of finding the boundary conditions between 
the different materials. Consequently, complex 
structure can be easily modeled simply by changing 
the property of elements In the calculatlonal 
domain. 

For a two media region consisting a part air 
and part bulk absorber, a model was used to describe 
the bulk absorber properties In two directions. An 
experiment to verify the numerical theory was con- 
ducted In a rectangular duct with no flow and with 
absorbing material mounted on one wall. Changes In 
the sound field, consisting of planar waves, was 
measured on the wall opposite the absorbing mate- 
rial. As a function of distance along the duct, 
fairly good agreement was found In the standing 
wave pattern upstream of the absorber and In the 
decay of pressure level opposite the absorber. 

The numerical formulation Is relatively simple 
to use and appears to give accurate modeling of 
experimental data. The theory may be a useful tool 
in the evaluation of viscous loss coefficients In 
bulk materials. 
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TABLE 1 


Position 

x/3.81 

Pb/PA. 

rms 

Position 

x/3.81 

Pb/PA. 

rms 

1 

0.00 

1 .0 

25 

6.417 

0.1793 

2 



26 

6.667 

.0990 

3 

.5 

1.105 

27 

6.917 

.0723 

4 

.75 

1.078 

28 

7.167 

.1304 

5 

1 .00 

.9904 

29 

7.417 

.1983 

6 

1.25 

.8655 

30 

7.667 

.2541 

7 

1.5 

.7066 

31 

7.917 

.3009 

8 

1.75 

.5791 

32 

8.167 

.3383 

9 

2.00 

.5358 

33 

8.833 

.3944 

10 

2.25 

.6013 

34 

9.083 

.3981 

11 

2.50 

.7428 

35 

9.333 

.3929 

12 

2.75 

.887 1 

36 

9.583 

.3790 

13 

3.00 

1.004 

37 

9.833 

.3582 

14 

3.25 

1 .077 

38 

10.083 

.3323 

15 

3.50 

1 .088 

39 

10.333 

.3141 

16 

3.75 

1.049 

40 

10.583 

.3078 

17 

4.417 

.7495 

41 

10.833 

.3155 

18 

4.667 

.6395 

42 

11.083 

.3332 

19 

4.917 

.5510 

43 

11.333 

.3563 

20 

5.167 

.4882 

44 

11.583 

.3817 

21 

5.417 

.4386 

45 

11.833 

.3934 

22 

5.667 

.3900 

46 

12.083 

.3960 

23 

5.917 

.3296 

47 

12.333 

.3894 

24 

6.167 

.2599 

48 

12.583 

.3738 
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(A) LOCALLY REACTING ABSORBING ADMITTANCE BOUNDARY 
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(B) EXTENDED REACTING INHOMOGENEOUS ABSORBING BOUNDARY. 

FIGURE 1. - BOUNDARY CONDITIONS EMPLOYED TO SIMULATED SOFT 
WALL WIND TUNNEL. 



THEORIES AT 
INTERFACES 

FIGURE 2. - TWO DIMENSIONAL DUCT FINITE ELEMENT MODEL. 
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FIGURE 3. - NO-FLOW ACOUSTIC DUCT TEST SECTION AND INSTRUMENTATION. 
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NORMALIZED RMS PRESSURE 



AXIAL CORRDINATE, Z 

FIGURE 5. - A COMPARISON OF THE MAGNITUDE OF THE AXIAL 
RMS PRESSURE VARIATION ALONG THE LOWER WALL IN A UNI- 
FORM DUCT WITH HARD WALLS AND A CHANGE IN PROPERTIES 
OF THE MEDIA AT Z OF 0.25 = 1.0 AND £ 2 = 4. 0) 

AS OBTAINED BY USING AN EXACT SOLUTION AND A FINITE 
SOLUTION FOR A PLANE WAVE (MODE-ONE) INCIDENT AT 
Z = 0 WITH U = 2tt. 
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AXIAL COORDINATE, X 


FIGURE 6. - EXPERIMENTAL AND THEORETICAL AXIAL PRESSURE 
PROFILE WITH PROPERTIES FROM HERSH's MODEL (U = 1.047). 
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FIGURE 7. - EXPERIMENTAL AND THEORETICAL AXIAL PRESSURE 
PROFILE WITH PROPERTIES FROM MODIFIED HERSH's MODEL 
(g = 1.047). 
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